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Abstract. We discuss the formation, evolution and observational param- 
eters of the population of short-period (^SlOhr) low-mass black-hole binaries 
(LMBHB). Their evolution is determined by the orbital angular momentum loss 
and/or nuclear evolution of the donors. All observed semidetached LMBHB are 
observed as soft X-ray transients (SXTs). The absence of observed short-period 
stable luminous X-ray sources with black holes and low-mass optical compo- 
nents suggests that upon RLOF by the donor, the angular-momentum losses 
are substantially reduced. The model with reduced angular-momentum loss rea- 
sonably well reproduces the masses and effective temperatures of the observed 
secondaries of SXTs. Theoretical mass-transfer rates in SXTs are consistent 
with those deduced from observations only if the accretion discs in LMBHB are 
truncated. The population of short-period LMBHB is formed mainly by systems 
which at RLOF had unevolved or slightly evolved donors (abundance of hydro- 
gen in the center X c ^0.35). Our models suggest that a very high efficiency of 
common envelopes ejection is necessary to form LMBHB. 



1. Introduction 

Currently, ten Galactic dynamically-confirmed black-hole candidate X-ray bina- 
ries with K/M spectral type secondaries and orbital periods ^1 day have been 



observed (jMcClintock Sz Remillardl 120061 ) . All these objects are SXTs. Their 



X-ray luminosity may vary by 5-8 orders of magnitude between quiescence and 
the peak of the outburst. Recurrence times spread from a about a year to tens or 
years and could be even longer. Their variability is interpreted as resulting from 
a thermal-viscous instabili ty of irradiat ed accretion discs around black holes in 
semidetached binaries (see Lasota 200 ll . and references therein). The estimated 



number of low-mass black-hole binar ies (LMBHB) in the Galaxy ra nges from 
several hundred to several thousand ( Chen et al.lll997l ; lRomanl [l998). The list 



of short-period black-hole SXTs and some of their observed and inferred param- 
eters are rendered in Table 1. Obse rvational data presented in Table 1 is based 
on the survey of the literature (see lYungelson Lasotal 120081 ) . 



Below, we consider the formation, evolution and some observational prop- 
erties of LMBHBs. 
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Table 1. Known short-period black- hole SXTs. 



No. 


Name 




Sp 


1 IK T \ 

(M recc ) 




M in , 








hr 




M Q yr" 


l 


M Q yr _1 




1. 


XTE J1118+480 (KV UMa) 


4.10 


K5-M1 


1.9 x 10" 


VI 


3.0 x 10~ 


10 


2. 


GRO J0422+32 (V518 Per) 


5.09 


M2-M4 


1.3 x 10" 


n 


1.9 x 10" 


10 


3. 


GRS 1009-45 (MM Vel) 


6.84 


K7-M0.5 


4.4 x 10" 


n 


1.8 x 10" 


10 


4. 


XTE J1650-500 


7.68 


K4 


2.0 x 10" 


n 


2.8 x 10" 


10 


5. 


A0620-00 (V616 Mon) 


7.75 


K3-K7 


3.3 x 10" 


n 


2.8 x 10" 


10 


6. 


GS 2000+25 (QZ Vul) 


8.28 


K3-K6 


2.0 x 10" 


10 


5.5 x 10" 


10 


7. 


XTE J1859+226 (V406 Vul) 


9.12 


G5-K0 


4.1 x 10" 


10 


1.0 x 10" 


9 


8. 


GRS 1124-68 (GU Mus) 


10.39 


K3-K7 


3.4 x 10" 


10 


3.5 x 10" 


10 


9. 


H 1705-25 (V2107 Oph) 


12.50 


K3-K7 


5.5 x 10- 


11 


4.1 x 10" 


10 


10. 


4U 1543-47 (IX Lup) 


27.0 


A2 


4.2 x 10" 


10 


1.3 x 10" 
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Note: (M recc ) - mass-transfer rate estimate based on recurrence times, M- m — upper 
limit to mass-transfer rate based on assumption of maximal truncation of accretion 
discs in SXTs. 



2. The model 

To obtain a model of the population of LMBHBs one needs to follow two steps: 
(i) the time-dependent formation of the population of detached "black-hole + 
main-sequence star" binaries and (ii) the subsequent evolution of every binary 
till the Hubble time. 

The threshold for the masses of black-hole prod ucing binary components 
is (20 - 25)Mq (e.g., Ergma van den Heuvellll998l ). The masses of the sec- 



ondaries in SXTs suggested by their spectral types are ^1.5M . Therefore 
the sequence of transformations of a binary which results in the formation of a 
LMBHB may be the following: 

(a) the primary component evolves off main-sequence and becomes a supergiant, 

(b) the supergiant overflows Roche lobe and forms a common envelope, since 
mass ratio of components q S> 1; if the components do not merge, the primary 
becomes a Wolf-Rayet (WR) star, 

(c) the WR star explodes as a supernova and forms a black hole. 

If the binary is not disrupted by the supernova explosion, a system with a 
black hole accompanied by a low-mass main-sequence star, i. e., a LMBHB, is 
formed. If the separation of components is sufficiently small (^1OR ), the orbital 
angular momentum loss (AML) via magnetic braking (MB) and/or gravitational 
wave radiation (GWR) may bring secondary component of the system to the 
Roche- lobe overflow (RLOF). 

The range of post-common-envelope binary separations which allow the 
formation of semidetached LMBHBs is very narrow. Hence the probability of 
formation of a SXT depends very strongly on the physical processes that de- 
termine the semi-major axis of the binary, i.e. the stellar winds at all stages of 
evolution and, especially, the efficiency of ejection of the common envelope. For 
the latter, the ratio of final aj to initial separations of components is equal 
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to (|WebbinkHl984l ) 



Of 



Mi, c 
Mi 



1 + 



a ce Ari r 



'Mi - Mi,, 
. M 2 



n -1 



(1) 



where a ce is the common envelope ejection efficiency, A - the parameter of the 
binding energy of the stellar envelope, Mi and Mi )C are the initial mass of the 
mass- losing star and the mass of its remnant, ri £, is the dimensionless radius 
of the star at the beginning of mass transfer, M2 is the mass of companion. If 
Mi 3> M2, then, crudely, df/at oc a re \- Both terms in the latter expression 
are h i ghly uncertain (see discuss ion in Podsiadlowski et al.ll2003l ; iKiel <fe Hurley 
20061 ; Yungelson Lasotal 120081 ) . Indirect estimates involving formation sce- 
narios for binaries with neutron-star or black- hole components suggest a ce \^2. 
This might mean that sources other than the orbital energy are involved in the 
ejection of the common envelope, though, the nature of these sources is still not 
fully understood. In practice, a ce X remains a parameter, tuning of which allows 
reproducing the properties of specific systems or stellar populations. 

In the second step of modeling, the evolution of every system is followed 
taking into account the AML via gravitational wave radiation and/or magnetic 
stellar wind (MSW) and, if necessary, nuclear evolution of the main-sequence 
star. 



3. The galactic population of LMBHBs 



In lYungelson et al.1 (j2006T ) and lYungelson fc Lasotal |2008) was carried out the 
population synthesis for galactic LMBHBs for three values of a ce X: 0.1, 0.5 and 
2. We refer the reader to these papers for the details of the computations. Two 
main conclusions were drawn from the models computed in these studies. 

First, the model for a ce X =0.1 is not compatible with observations of SXTs 
since in this case the overwhelming majority of black holes has masses ^,14 M , 
exceeding the l argest estimated bl ack hole mass in known SXTs (9.7 ± 0.6 for 



Froning et~aLll2007l ). 



A0620-500, see 

Second, if the AML via MB for LMBHBs is treated in a 



"standard" way, 

assuming after Verbunt Zwaanl ( 1981) that for com ponents of close binaries 
the braking law for single field stars ( Skumanich 19721 ) can be extrapolated over 
the range of rotational velocities from several 10 to several 100 km/s and that 
the spin-orbit coupling is efficient, model mass-transfer rates for LMBHBs at 
-Porb~2hr are so high that these systems might have stable hot disc s according 
to the disc instability model (DIM) criterion of iDubus et al.l ( 19991 ). But such 



stable and bright LMBHBs have not been observed. 

A strong reduction of the magnetic braking efficiency in close binaries as 
compared to s ingle stars is suggested also by the data on ste llar rotation in young 
open clusters ( Collier Cameronl 20021 ; lAndronov et al.ll2003l) and the mass trans- 



fer ra tes in cataclysmic variable s (e.g ., Hameurv et al.Tl988l; llyanova fe Taaml 



2003). In lYungelson et al.l (|2006l ) and lYungelson fc Lasotal ^20081 ) we computed 



a population model of LMBHBs assuming that MB stops operating once the 
RLOF occurs ("no-MB" model). We found that the number of galactic LMB- 
HBs remains of the same order of magnitude as in the case with active MB 
(several 1000) but all systems are transient. 



4 L.R. Yungelson and J.-P. Lasota 




5 10 15 5 10 15 

P, hr P, hr 



Figure 1. Left panel: Model population (dots) vs. observational estimates 
of the ranges of effective temperatures of donors in SXTs inferred from their 
spectral types (vertical lines). The latter are annotated according to their 
number in Table [T] Large filled circles give T c ff of donors derived from the 
fits to synthetic spectra. Heavy solid lines to the left and right, respectively, 
show the approximate limits of the region which may be occupied by LMB- 
HBs, while thin solid lines show tracks for the donors in (1.1 + 4)M Q systems 
which have at RLOF X c =0.45 and 0.37, respectively (see the text for discus- 
sion) . The Sp — T c ff relation used in the paper is shown at the right border 
of the coordinate box. Right panel: Masses of donor-stars in modeled popu- 
lation. Vertical lines show the ranges of M2 inferred from their spectral types 
(annotation like in the left panel). Heavy and thin solid lines - the same 
tracks as in the left panel. 



In lYungelson et all ( 20061 ) and lYungelson &: Lasotal ( 20081 ) we compared the 
a ce \ =2 "no-MB" model with observations. Below, we present the model for 
a ce X =0.5. In this case, the number of LMBHBs that reach contact in Hubble 
time and evolve to -P or b~10hr is ~ 6000 and ~ 3000 of them have currently 
P orb >1.5hrE| For comparison, in the model for a ce X = 2 these numbers are 
about 12000 and 5000, respectively. 



3.1. Effective temperatures and masses of secondaries in LMBHBs 

When spectral types of the secondaries in observed SXTs are known, one can 
evaluate their effective temperatures and masses and to compare them with 
model predictions. Inevitably, such a comparison may be only crude, since the 
spectra of the secondaries are contaminated by the radiation of accretion disks 
and hot spots and we have at our dispo sal Sp — T P ff and Sp — T e g relations for 
main-sequence stars only (relations from lCoxll200Cll . were used). In the left panel 
of Fig. [T] we plot the effective temperatures of donor-stars in modeled systems as 
a function of orbital period and compare them with the ranges of T e g inferred 



1 Systems with shorter P or b have and the character of m ass exchange in them is unclear 

because of the effect of resonances, see Yungelso n et alj (|2006T ) for details. 
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from the spectral types of donors in SXTs. Having in mind uncertainties in 
the spectral type determinations and Sp — T e g- scale, we find that the model 
satisfactorily reproduces T e g of the donors in the LMBHBs with P rb~9 hr. 

In our models the efficiency of magnetic braking is 0. In reality some MB 
might b e acting as, probably, i n the case of cataclysmic variables (see, for in- 
stance, llvanova <fe Taaml 120031 . but also 13.2.1 below). We show in Fig. Q] two 
"limiting " tracks: for the system with initial masses of components 1 and 12 
Mq and post-circularization period -P O rb,0 = 0-4 day in which the donor is al- 
most unevolved at RLOF and MB does not act after RLOF and for a (1+4) Mq, 
-Porb,o = 1-9 day system in which at RLOF donor has hydrogen abundance in 
the center X c ~ 1CP 4 and the MB continues to operate dur ing mass-transfer 
with an efficiency corresponding to IVerbunt &: Zwaanl ( 1 98 ll ) law. Crudely, if 
the efficiency of MB is not 0, the model population must be located between 
these two limiting curves. Of course there will be a contribution from lower and 
higher mass systems, as we plotted the 1 Mq tracks for simplicity only. Adding 
some MB to our model will shift the population to the right, providing a better 
agreement with observations. Adding moderate MB will influence mainly the 
long-period systems, without producing stable systems. 

A similar satisfactory agreement of model populations with observations is 
found for the masses of secondaries of SXTs (Fig. [Q right panel) . 

We note that, since for large initial q the transformation of separation of 
components in common envelopes depends linearly on a ce X, the distributions of 
model populations in P — T e jf and P — M2 plots are similar for a ce A=0.5 and 
2, and only differ by the density of points per unit area of the diagrams. 

There are SXTs with P Qrb « (8-12) hr - GRS 1009-45, XTE 1650-500, 
A0620-00, and GS 2000+25, which, apparently, are located below the "popu- 
lated" area in Fig. [TJ However, we should note that in our modeling we tried 
to avoid the effect of bifurcation of evolutionary tracks - evolution of systems 
to shorter or longer periods upon RLOF, depending on the extent of hydrogen 
depletion in the cores of the models. For this reason we restricted ourselves to 
systems evolving to shorter P or b only. However, our computations show that, if 
GWR is the sole sink of angular momentum, binaries with M\ ~ (4—12) Mq and 
M<it^A Mq evolve to longer periods if at the instant of RLOF their secondaries 
have X c ^,0A. This is illustrated in Fig. Q] by evolutionary tracks for initially 
(1.1+4) Mq systems in which secondaries have, respectively, X c ~ 0.45 and 0.37 
at the RLOF. The latter binary spends several Gyr evolving to longer periods. 
Initial (post-circularization of the binary after supernova explosion) orbital pe- 
riods of these two binaries Po differ by 0.1 day only. Since the distribution of the 
binaries over Pq is continuous, this proximity of initial parameters of the bina- 
ries and striking difference in their evolutionary behavior suggests a possibility 
of explaining the origin and parameters of SXTs with P or b ~ (8 — 12) hr. This 
conjecture has yet to be confirmed by detailed modeling. 



3.2. Mass-transfer rates in LMBHBs 

An estimate of mass-transfer rate in an SXTs may be obtained by dividing the 
mass accreted during outburst by the recurrence time. However, (i) recurrence 
times are known for A0620-00 (about 60 yr) and 4U 1543-47 (about 10 yr) only; 
(ii) it is not evident that the rate calculated this way represents the secular 
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Figure 2. Mass-transfer rates in model LMBHBs as a function of their 
orbital periods (dots). Arrows mark upper limits to the estimates of mass- 
transfer rates in observed SXTs as given by Eq. © . Crosses are estimates of 
mass-transfer rates in SXTs based on recurrence times. 



value, and (iii) this method of estimate assumes that between outbursts, when 
accretion disc is "refilled", accretion onto black hole do es not occur. Th e last 
assumption is put in doubt both by observations (see e.g. iDone et alJl2007l . and 
references therein) and models (see lLasotall2008l . and references therein) which 
suggest that quiescent discs in SXTs are truncated and therefore leaky. In the 
latter case, mass-transfer rate cannot be l arger than the critica l-for-stability 



accretion rate at the truncation radius (see lYungelson et al.l l2006. for details). 



The actual mass-transfer rate should be between the values estimated by the 

two methods. 

(120081 ): 



Using the expression for the critical accretion rate lLasota et al 



2.64 • 10 15 a° f i?? n 58 M 



-0.85 



40 



'1 



(2) 



where a is viscosity parameter and R\o is disc radius in units of 10 10 cm, we 
obtain for the upper limit of mass-transfer rate 



M<2.5 ■ W~ 7 



(0.5 -0.227 log . 



10.32 



Py 2 fF° 8 M yr-\ (3) 



where P^ - orbital period in days, /i^,0.48 - fractional disc truncation radius. 

In Fig. [2] we compare the model mass-transfer rates with the estimates 
obtained by the two above mentioned methods. The "leaky disc" estimates of 
M max for XTE J1118+480, GRO J0422+32, and GRS 1009-45 strongly suggest 
that the AML in short-period LMBHBs might be really defined by GWR only. 



4. Conclusion 

We calculated models of the Galactic population of short-period low-mass black- 
hole binaries which are identified with soft X-ray transients. We found that 
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using the values of the common-envelope parameter ot ce \ between ~ (0.5 — 2) 
and assuming a strongly reduced magnetic braking it is possible to reproduce, 
(within the uncertainty of observations) the number of the LMBHBs in the 
Galaxy and the effective temperatures and masses of the donors in these systems 
(as inferred from the spectra of the stars). The above mentioned values of a ce X 
imply that the common-envelope expulsion in the progenitors of SXTs has to 
be very efficient and that sources other than orbital energy may be required in 
this process. 

In our model, all short-period LMBHB systems are transient in agreement 
with observations. 

Model mass-transfer rates in LMBHBs are consistent with the upper limits 
derived from observations under assumption that accretion discs in SXTs are 
leaky. 
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